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The domain patterns in a thin ferromagnetic film are investigated in both experiments and 
numerical simulations. Magnetic domain patterns under a zero field are usually observed after 
an external magnetic field is removed. It is demonstrated that the characteristics of the domain 
patterns depend on the decreasing rate of the external field, although it can also depend on other 
factors. Our numerical simulations and experiments show the following properties of domain 
patterns: a sea-island structure appears when the field decreases rapidly from the saturating 
field to the zero field, while a labyrinth structure is observed for a slowly decreasing field. The 
mechanism of the dependence on the field sweeping rate is discussed in terms of the concepts 
of crystallization. 
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There are a large number of physical and chemical 
systems that display domain patterns. 1 Ferromagnetic 
garnet thin films also show domain patterns, which can 
have various kinds of structures. For example, when an 
external field is applied, a hexagonal lattice, bubbles, 
or stripes are observed. Under a zero field, usually a 
labyrinth structure is observed, although a parallel-stripe 
structure is more stable than a labyrinth structure. In 
fact, it is observed in experiments that a labyrinth struc- 
ture changes to a parallel-stripe structure under field cy- 
cles. ' However, the domain patterns observed under 
a zero field do not always have a labyrinth or parallel- 
stripe structure. In some samples, under a certain con- 
dition, a sea-island structure is observed, which consists 
of many small up-spin (down-spin) domains surrounded 
by a "sea" of down (up) spins. 

In order to observe the magnetic domain patterns un- 
der a zero field in experiments, usually an external mag- 
netic field is applied once up to the saturation field, and 
then it is removed. Here, we will focus on the domain 
patterns under a zero field and their dependence on the 
decreasing rate of the external field, although it is nat- 
ural that the characteristics of the observed patterns 
can depend on the properties of the samples, e.g., dis- 
order, anisotropy, interactions among spins, and so on. 
In fact, there have been several simulations that show 
how domain patterns change depending on such prop- 
erties under a very slowly changing field. 4-6 There also 
exist some works on the dependence of hysteresis on the 
frequency and amplitude of an oscillating field. 7-9 How- 
ever, for the characteristics of domain patterns under a 
zero field, their dependence on the changing rate of a 
time-dependent field has not been investigated so far. 

In this letter, we report both the experimental observa- 
tions and the numerical simulations of domain patterns. 
In particular, we show how the characteristics of domain 
patterns depend on the sweeping rate of the magnetic 
field. First, the experimental procedure and results are 
explained. Next, the model of the system and numeri- 
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Fig. 1. Experimentally observed domain patterns under a zero 
field. The external field was decreased at a (a) fast rate and (b) 
slow rate. (See the text.) 



cal procedure are described. After displaying numerically 
simulated domain patterns, we study the field depen- 
dence of the characteristics of domain patterns, i.e., the 
number of domains and the domain area. The mechanism 
of the dependence will also be discussed. 

Experiments were performed on ferrimagnetic garnet 
films of composition (BiGdY) 3 (FeGa)sOi2 grown epi- 
taxially on the (111) face of gadolinium gallium gar- 
net substrates. This film has strong uniaxial magnetic 
anisotropy with the easy axis perpendicular to the film. 
The perpendicular saturating field is about 120 Oe. Mag- 
netic domain structures are easily observed by using a po- 
larizing microscope with the help of the magneto-optical 
Faraday effect. Figure 1(a) shows the domain pattern 
under a zero field when the magnetic field has been de- 
creased from 220 to Oe at a fast sweeping rate of 2 x 10 5 
Oe/s. This picture shows a sea-island structure. The ob- 
served structure at a sweeping rate of 10 Oe/s is shown 
in Fig. 1(b). At this slow rate, the domain pattern has a 
labyrinth structure. 

In our numerical simulations, we apply a simple two- 
dimensional Ising-like model (see refs. 4, 5, and refer- 
ences therein). We consider a scalar spin field <fi(r), where 
r = (x, y). Namely, the model is a continuous model. The 
positive and negative values of 4>{r) correspond to the 
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up and down spins, respectively. Our model Hamiltonian 
consists of four energy terms: uniaxial-anisotropy energy 
iJ an i, local ferromagnetic interactions (exchange interac- 
tions) Hj, long range dipolar interactions and in- 
teractions with the external field H cx . The anisotropy 
energy term implies that 4>(r) tends to take either one of 
two different values [<j>(r) = ±1]: 



a\(r) 



:lr | ^) 2 , <t>(rf 



(1) 



where the factor A(r) expresses disorder, and 

A(r) = l + /i(r)/4. (2) 

Here, /u(r) is an uncorrelated random number with a 
Gaussian distribution whose average and variance are 
and /1q, respectively. However, /u(r) should have a cutoff 
so that A(r) is always positive. We consider the disorder 
effect only in the anisotropy term, although some effects 
of disorder could also be applied to the other interaction 
terms. The ferromagnetic (exchange) and dipolar inter- 
actions are described by 

|V0(r)| 2 



Hj = dr 



and 



ff di =7 J drdr'((>{r)(i>{r')G{r,r'), 



(3) 



(4) 



respectively. Here, G(r, r') ~ \r — r'\~ 3 at long distances. 
The integral in eq. (4) should have a lower cutoff since 
G(r,r') — > oo at short distances. The term from the 
interactions with the external field is given by 



-h(t) J <W(r), 



(5) 



where h(t) is the time-dependent external field. From 
cqs. (1) (5), the dynamical equation of our model is de- 
scribed by 

d(j>{r) 5{H &ai + Hj + H di + H cx ) 



dt 



5<p(r) 

a\(r)[<P(r) - </>(r) 3 ] + 0V 2 <p(r) 

J dr'<t>(r r )G(r,r') + h(t). (6) 
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The coefficients of the anisotropy energy and the ex- 
change and dipolar interactions — a, 0, and 7 respec- 
tively — are all positive. 

In numerical simulations, it is useful to calculate the 
time evolution of eq. (6) in Fourier space. The equation 
is rewritten as 

^ = a[X{</> ~ <f)\k ~ (0k 2 + 7 G fc )0 fc + h(t)S kfi , (7) 

where [•]& denotes the convolution sum and Gk is the 
Fourier transform of G(r, 0). In the calculation below, if 
G(r,0) eee 1/H 3 , then 



Gfe = a - aifc, 



where k — \k\ and 



ciq = 2ir I rdrG(r), a\ 
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Here, d is the cutoff length, which is approximately the 
thickness of the film in experiments. We set d = 7r/2, 
namely, ao = 4, in our simulations. In the simulations, 
512 x 512 lattices with periodic boundary conditions 
are used, and we fix two parameter values: = 2.0, 
7 = 20/ai = 2.0/71". These parameters determine the 
characteristic length of domain patterns. According to 
the linear analysis of eq. (7), <pk of the modes with 
k = 1.0 grows the best at these fixed values. Namely, 
the characteristic wave number fco of the domain pat- 
terns is unity: fco = 1.0. Moreover, we fix the other two 
parameters as follows: a — 2.5 and [i = 0.3. It would 
be interesting to investigate the dependence on a or /io, 
which will be studied elsewhere. Only the changing rate 
v of h(t) is the control parameter of our interest: 



h(t) = h ini - vt, 



(10) 



where /i; n i is the saturating field; h lnl ~ 1.4 in this case. 
The value /i; n i was estimated from a simulation under an 
ascending field. In other words, the domains first spon- 
taneously form a pattern under a zero field, and then an 
increasing field is applied. At a certain value of h, the do- 
mains with negative (f>(r) disappear, and we define this 
value as /i; n i- 
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Figure 2 shows the domain patterns under a zero field 
that are simulated by the abovementioned method. The 
white and black areas correspond to positive and neg- 
ative </>(r), i.e. up and down spins, respectively. As the 
initial condition, <f>(r) was given by random numbers in 
the interval 1.0 < 4>{r) < 1.1 at t — 0. The external field 
decreases from ft,; n ; to 0, following eq. (10). Once the ex- 
ternal field becomes zero at i = h- m \/v, the field remains 
zero for t > t . The change in the domain patterns for 
t > t is very little in this case. For some other sets of 
parameters (a and /in), however, the characteristics of 
the domain patterns change after the field vanishes. 10 

The simulation results have the same behavior as the 
experimental ones. Namely, a sea-island structure, which 
consists of many short domains, appears under a rapidly 
decreasing field, and a relatively small number of do- 
mains grow to form a labyrinth structure under a slowly 
decreasing field. In other words, both the experimental 
and numerical results suggest that the field sweeping rate 
affects the domain size (domain length). 

The abovementioned results also mean that a domain 
pattern forms an unstable structure under a rapidly de- 
creasing field. Since the ferromagnetic (exchange) inter- 
actions prefer small domain surfaces (domain walls), a 
structure composed of a few large domains is more sta- 
ble than one composed of many small domains. In other 
words, a labyrinth structure is more stable than a sea- 
island structure from the viewpoint of energy. Before dis- 
cussing why such an unstable structure appears under a 
rapidly decreasing field, we consider how domains grow 
to form a structure under a descending field. 

How are the domain patterns in Fig. 2 created? The 
field dependence, which is now identified with the time 
dependence due to eq. (10), of the number and area of 
negative-^ (black) domains is shown in Fig. 3. Under a 
descending field, the total number of domains for each v 
rapidly increases at a certain field and has a plateau in 
the low-field region. In the plateau region, each domain 
grows; however, no new domain appears. The length of 
the plateau is longer for small v than for large v. For the 
rapidly decreasing field (v = 10~ 2 ), the value of the field 
where the first domain appears is lower than that for the 
other slowly decreasing fields (y = 10 -3 , 10~ 4 ). From 
these results, we see that a small number of domains 
grow for a long time under a slowly decreasing field and 
vice versa. 

Both the average and total domain areas show rapid 
growth at first and slow growth in the low-field region. 
This behavior is explained as follows. In the rapid growth 
region, where the domain area increases rapidly, domains 
grow freely. In the slow growth region, however, domains 
cannot grow freely since they tend to repel each other be- 
cause of the dipolar interactions. Moreover, domains do 
not connect with each other for the parameter combina- 
tion used in this calculation. 11 In the small-i> case, since 
the number of domains is small, the change in the growth 
rates becomes clear. Here, we note that the difference 
between the total domain area and 512 2 /2 in Fig. 3(c) 
corresponds to the magnetization. We see that there is a 
small difference between the total domain area under a 
zero field and 512 2 /2, which corresponds to the remanent 
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Fig. 3. The field dependence of (a) total number of domains, 
(b) average domain area, and (c) total domain area. The solid, 
chained, and dashed lines are for v = 10~ 2 , v = 1Q~ 3 , and 
v = 10~ 4 , respectively. The domain area is the number of grid 
points that have negative values of <j>(r). The results under a zero 
field correspond to Fig. 2. 



magnetization. 

Now, let us discuss why the characteristics of domain 
patterns depend on the field sweeping rate, as shown in 
Figs. 2 and 3. We can explain the mechanism of this 
behavior by using the concept of crystallization: The 
driving force which causes nucleation corresponds to su- 
persaturation or supercooling, and the nucleation energy 
AF is small for high supersaturation (supercooling). 12 In 
this system, the nucleation energy AF is the energy that 
is needed to create a negative- </> (black) domain in the 
"sea" of positive-^ spins (white space). However, there is 
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no denned quantity called supersaturation or supercool- 
ing. Let us consider a quantity that expresses the degree 
of non-equilibrium as "supersaturation." The supersatu- 
ration should depend on the external field and the field 
sweeping rate. Since the changing rate of the valuable 
4>(r) is finite, the system cannot always keep up with the 
change in the field, ft is natural to expect the following: 
the larger the value of v, the higher is the supersatura- 
tion. Therefore, when v is large, AF is small because of 
high supersaturation and nucleation easily occurs, which 
leads to the sudden appearance of many nuclei. On the 
contrary, when v is small, only a small number of nuclei 
appear because of large AF. Once the negative-^ do- 
mains appear, the domains grow but new domains hardly 
appear. 

The supersaturation as well as the nucleation energy 
AF should be related to the hysteresis loop. In rcf. 9, it 
is shown that the shape of the hysteresis loop changes 
depending on the sweeping rate of the field. The area 
of the hysteresis loop is equal to the energy dissipated 
per cycle. The area is large for a fast sweep and small 
for a slow sweep. 9 Though the hysteresis and dynamic 
phase transition were mainly studied for different models 
in refs. 7-9, the methods used in the references will be 
helpful to evaluate AF. The quantization of AF will be 
discussed sometime in the future. 

There is another question to be considered: how is it 
that a labyrinth structure is observed in some samples 
even when the external field is suddenly removed? This 
can be explained by the difference in the proper time 
scale of the individual samples in the absence of the ap- 
plied field. The proper time scale is determined by the 
anisotropy, interactions between spins, and disorder in 
the samples. The difference in the time scale may be 
interpreted as the difference in "the sensitivity to the 
dynamical change in the field." Using the concept of the 
sensitivity to the field change, we can consistently under- 
stand both our results and the above fact that a labyrinth 
structure often appears under a sudden decay of the field. 
Our results suggest that the sensitivity, in the same sam- 
ple, should depend on the field sweeping rate v. On the 
other hand, if a sample has very high sensitivity to the 
field change, a labyrinth structure can be observed in the 
sample even when the external field is suddenly removed. 
In other words, the balance between the sensitivity and 



the field sweeping rate is one of the factors that causes a 
variety of domain patterns. 

In conclusion, we have confirmed that the characteris- 
tics of the magnetic domain patterns under a zero field 
depend on the decreasing rate of the external field in 
both the experiments and the numerical simulations. For 
a rapidly decreasing field, many short domains form a 
sea-island structure. On the contrary, a small number 
of long domains form a complex labyrinth structure for 
a slowly decreasing field. Those properties are explained 
qualitatively by using the idea of crystallization. In other 
words, in the case of large v, high "supersaturation" low- 
ers the nucleation energy, which causes a sudden appear- 
ance of many nuclei, i.e., negative-;/) domains, in a sea of 
positive-0 spins. However, a labyrinth structure appears 
even when the external field is suddenly removed in some 
experiments because the spins in actual garnet films usu- 
ally vary fast enough to keep up with the change in the 
field. The results in this letter show that the characteris- 
tics of domain patterns are dependent on the field sweep- 
ing rate, or more generally, on the balance between the 
rate and a dynamical property such as "the sensitivity 
to the dynamical change in the field." 
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